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I . INTRODUCTION 



The study of the dynamic stability of airplanes "by 
the method of snail motions involves the determination of 
oscillatory and aperiodic motions. There are a great 
number of well-established methods that permit us to cal- 
culate the periods and the damping factors of these mo- 
tions but which do not, in general, give any indication of 
the amplitudes of the motions set up. We shall in the 
present paper, with the aid of a number of numerical ex- 
amples, attempt to clarify the phenomena arising after a 
series of typical initial disturbances. It will be found 
that the results of these calculations contribute to an 
understanding of the mechanism of the motions which are 
set up, and it is our belief that the publication of these 
results, although limited to particular cases, will facil- 
itate the study of the motions and flight paths. In what 
follows, we shall strictly separate the study of the lon- 
gitudinal from that of the lateral motions and shall not 
consider the complete theory. 



NOTATION 



The OX, OY, and OZ axes will be assumed fixed on the 
airplane (fig. l). The projections of the velocity on the 
axes will be denoted by u, v, and w, respectively. They 
define the angles of attack i, and of yaw j, these two 
angles fixing the position of the airplane on its path. 
When the components v and w are small with respect to 
u , we may write : 



w 

1 = - — 

u 



= + 2 

U 



"Etude theorique de quelques trajectoires d'avions apres 
perturbation initiale." Bulletin du Service Tech- 
nique de L 'Aeronautique, No* 17, June 1937. 
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and have, moreover, the relation 

V = y~u 2 + v 2 + w 2 a u 

The projections of the instantaneous rotational veloc- 
ity Q on the three axes are p, q, r, constituting the 
rolling, pitching, and yawing angular velocities. The 
rolling is positive if it tends to lower the right wing 
and raise the left; the pitching is positive if it tends 
to nose the airplane down, while the yawing is positive if 
it tends to turn the airplane toward the left. The con- 
trols corresponding to these rotations are, respectively, 
the ailerons, elevator, and rudder, their deflections a, 
3, V heing positive when they tend to turn the airplane 
in the positive sense of the rotations. 

The position of the system of axes OXYZ fixed to 
the airplane, will he determined as a function of a system 
of spatially fixed axes OX Q Y 0 Z Q hy means of the three 

angular, displacements defined in figure 2. The deriva- 
tives of these angles 

dcp i 
— — s cp 
dt 

4fi = e* 

dt 
dt 



are connected with the angular 
geometric relations 

P = cp ' cos G - \|/' 

q = 6 ' + \|/ ' 

r = cp ' sin 8 + \|/ ' 



velocities p, q, r "by the 

cos cp sin 8 
sin cp 

cos cp cos 6 
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II. LONGITUDINAL HOTION 
Method Used 



The equations of motion referred to axes fixed to the 
airplane are, for nonuniform flight (fig. 3)* 



T + X + P sin 8 = ~ (~ + qyr} 

g \dt ' J 

z ^ p cos e ^) 

fir \dt / 



K 4- Ts 



g ^dt 

g ft* 

to which should he added the relation 



d9 
dt 



The equation £5 may he T7ritten: 

f l (u, W , q , 6 ) 



du 
dt 



dv? 
dt 



77 - ifl 



dt 



- f , 



d8 
dt 



(u, w, q , 6 ) 
(u, W, q, G ) 
(u, w, q, 6 ) 



r 



(I) 



(II) 



At the instant t 0 , let the values satisfying the 
system be u 0 , w 0 , q 0 , 8 0 , and let us give the airplane 
an initial disturbance, so that at the instant t 0 + 6t 
the variables have the values u 0 + £u, vr 0 + 8t7, etc. 
We nay replace the four variables by their increments Su, 
8w, 8q, 56 with respect to their initial values. This 



* In the section on longitudinal motion, r denotes the 
radius of gyration with respect to OY. 
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substitution, if the increnents are assuaed to he snail, 
makes it possible to linearize the system of differential 
equations. We thus obtain: 



cl(6u) Bf x 9fj. df 1 di\ 

— - 8u + -- L 6w + ~- 5q + — - 6 6 



dt 

d(8w) 
dt 

d(8q) 
dt 

a(56_) 

dt"" 



d u 



3f 2 

■— - 6u + 



d w 
3f 



3q 
3f. 



36 



3f. 



3u 



3f 3 

— - 8u + 
3u 3 w 



3tt 
of 3 



- 8w + --a- Sq + -- 2 - 6 6 



36 



3f 3f - 

8w + Sq + 6 6 

3q 36 



> (HI) 



3f 



9f : 



— *• 6u + — 6w + -— * 6a + ~-~ 8 6 
3u 3t7 3q 36 



In this new system the disturbances 8u, 8w, Sq, 6 6 
are the variables and the partial derivatives are con- 
stants whose values are determined by the initial condi- 
tions. This system is linear and readily integrated, 
yielding the variables Su, Su, Sq, 86 as functions o,f 
the time after an initial disturbance. The integrated 
system is of the form: 



8u 
8 w 

6q 
86 



C 4 e X 4 t 



x.t 



m 1 C 1 e A i t + ffi a 0 s e X 3* + n 3 C 3 e X 3 t + n 4 C 4 e X 4 t 
a x %e Xlt + n a C 2 e A 2 * 4- n 3 C 3 e X 3 t + a 4 0 4 ©*** 



MIY) 



cliara 


system the four value 
cteristic equation of 


s of 

the a; 






3i\ 




3f x 


3 u 




3q 




36 


3f 2 
3u 


* - x 

0 V/ 


3f 2 
3q 




3f 2 
36 


3u 


3f 3 
3w 


3f 3 
3 q 


A 


3f 3 

36 


3f 4 


3f 4 


3f 4 






3u 




a~q 




36 



■ 0 



(v) 



1I.A.C.A. Technical Lie no ran dun Ho, 867 



The groups of coefficient s li&i'&ii ^2 n 2 n 2* etc ', 
still "be obtained in accordance with the nethod of Lagrange, 
hy means of any three of the four equations of system (VI), 
into which are successively introduced the four values of 
X : 

__L - X + I -A + n * n £h = 0 



9 u V 3w / dq SG 







n 


a. 










- x) + 


n 


aq 




+ 










n 


^3 q 






n 




3 T7 


+ 





af 3 , of, /af 3 a af 3 

_-5- + X — 2. + n — i - A I + n 

au aw saq / ae 

af 4 , af 4 df 4 ^af. 

_i. .l i _ — + a — - + n 



(VI) 

■ 0 



(*L - x) . 



B u 3 \7 3 q \se 

and the constants of integration $% , 0 2 , C 3 , and C 4 , the 
latter depending on the initial conditions chosen. 

The characteristic equation (V) Is generally written: 

X 4 + A x a 3 + A 2 X 2 + A 3 X + A 4 = 0 (VII) 

This equation, in general, possesses imaginary roots. If 
the four roots A are imaginary 

X = a ± oi 
equations (IV) may ho written: 



SU m 




sin 


(bt 


+ 




) 


+ 


ea't Pi t 


sin 


(b 't 


+ 


(Pi 


«) 


8w = 


e at p 2 


sin 


(bt 


+ 




) 


+ 


e a,t p 2 ' 


sin 


(b 't 


+ 




') 


6q = 


e at p 3 


sin 


(ot 


* 


^3 


) 


+ 


e a,t p 3 ' 


sin 


(b't 


+ 


^3 


0 


66 = 


e at p 

1 4 


sin 


("bt 


+ 


^4 


) 


+ 


e a't p i 

4 


sin 


(b »t 


+ 


^4 





(Till) 



In the above form they permit us easily to follow the 
ferent disturbances which the airplane undergoes. 



dif- 
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When the four roots of the characteristic equation 
are imaginary, the notion of the airplane consists of the 
supcrpor.it ion of two oscillatory motions. In each of 
those p represents the maximum possible amplitude, 
sin ("bt + cp) determines the sinusoidal form, cp being 
the phase shift, and b a value such that the period T = 
2*nr/b, e a (equals 1 for t » 0) represents the damping* 

Tie then proceed as follows: 

1. Knowing the characteristics of the airplane stud- 
ied, tre deternine its characteristic equation (V) and ob- 
tain the roots A. From these roots we can find the period 
I and the damping factor e at . These characteristics de- 
pend on the airplane and not on the initial disturbance 
considered. 

Z% UTe next deternine by equations (VI) the groups of 
factors I, n, n, these factors being similarly independ- 
ent of the initial disturbance considered. When the roots 
\ are conjugate imaginary the factors I, a, n will also 
be conjugate imaginary and we nay write: 

lj = a + Pi n x » Y * 81 m « € + cpi 

l 2 = a - pi m 2 « Y - gi n 2 = € - cpi 

3« now introduce the initial conditions of the no- 

tion, that is, the disturbance considered. The latter is 
expressee 1 by a real value of one or several of the variable 
8U, Sw , 8 q , 86, and from it we nay determine the constants 
of integration, which are conjugate imaginaries: 

0 ! - A + 3 i C 3 = A ■ + 3 'i 

C 2 - A - 3i C 4 = A< - B»i 

The quantity G 1 e^i t s G^e^a* will assume the follow- 
ing forms : 

0 x e^i t -f© a e^a' t = C x e f* 4 ' 1 * 1 )* + Ce^-' 01 ^ 



e 



at 



[C 1 (cos bt+i sin bt)fC 2 (cos bt-i sin bt 



- e at [(Gi + Cs) cos bt + i(C a -C a ) sin bt] 
m e at [?. cos A bt - 23 sin bt] 

- e 3 -^ [sin cos bt + so s cp x sin bt ] 



s e^ 5 



p x sin (bt + g> ) 
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with 



sm qp x 



= 2 ,/a 2 + B g 



2A 

P, 



- 2B 

cos co ~ —~ — 

p x 

Similarly, the part of the root l 1 0 1 e x i t + l 2 C 3 o X2t 
will assume the following forns: 

e at + l 8 0 a ) cos bt + i (\ l 0 1 - l a C 2 ) sin bt] 

e at [(2A a - 2B 3) cos bt - (2A 3 + 2B a) sin bt] 
e at o sin (bt + cd„ ) 



with 



2 

sin cp 



(Aa - Bg) 2 + (A3 + Ba)' 



- 2 (Aa - B8) 



2 



P 2 



cos cp. 



- 2 (A3 + Ba) 



ffe nay note that p may "be written: 

2 



2 V 



Tie see that the values of the amplitudes p and of 
the phase displacements cp depend on the initial condi- 
tions. For the same oscillation the ratios P 2 / P x • P%tP%* 
p 4 /p are, however, independent of the initial distur "banco 
and dcDcnd only on the characteristics of the airplane. 
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AIRPLANES STUDIED 



17 e shall study a series of airplanes differing only 
in the coefficient of static stability 8C M /3i. If i 

may he taken equal to - ~, the derivatives with respect 

to w figuring in the calculations are merely the deriva- 
tives with respect to the angle of attack i multiplied 

- i« We have shown elsewhere that the longi tudinal 

motion of an airplane is sufficiently determined "by : five 
geometric characteristics; seven aerodynamic characteris- 
tics; a characteristic depending on the engine propeller 
unit; and a characteristic depending on the weight. In 
our example, we have assumed the following values for 
the so clement s : 

1. G-eometric characteristics: 

I tt 2.3 meters 

r » 1,3 n 

s = 0. 

I ■ = 2.6 I 

S » = i S 
7 

M and S* "being the lever arm and the area of the horizon- 
tail surface, respectively. 

2 . Aerodynamic character! sties : 

The quantities C x and C 2 are the coefficients of 

the aerodynamic forces in the directions of the 
axes fixed to the wing. 

The angle of attack i of the state considered is 0 
C z = 0.40. 

Minimum drag-lift ratio at this angle of attack: 
P = 0.125. 

§?X s 0.005; ||i „ 0.065. 
di di 
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The damping of the pitching oscillations will he char 
actcrized by a coefficient which we shall write in the 

ac* 

form TJi as if the entire damping were produced by the 

action of the horizontal tail surface. In the above de- 
pression C' 2 would he the lift coefficient of the tail, 

i ! the real angle of attack of the tail. 

As a matter of fact, we shall take a numerical value 
somewhat larger than that corresponding to the derivative 
of the lift of the tail in order to take account of the 
damping introduced "by the parts of the airplane other than 
the horizontal tail surface. We shall take: 



dC « 

z = 0.065 

di ' 



Finally, we shall take 



di 



the coefficient of static stability, as the variable ele- 
ment in the example. 

As the variations of the static stability mry be ob- 
tained by a simple displacement of the center of gravity 
without varying the tail surface, we may assume it to be 
possible to vary \i without varying the other character- 
istics. 

3 # Element characterizing the engine-propeller unit* 
Wc write : 

AX . ^ h 4 

dV V 
giving, h the value 0.5. 



4. Element characterizing the weight of the airplane 

The weight per square meter will be taken as equal to 
40 kilograms. From these data the result follows immedi- 
ately that the speed of the airplane under the conditions 
considered is 40 meters per second (90 miles per hour). 
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The above elements permit us to determine the charac- 
teristic "biquadratic equation in ft*, 

CALCULATED OSCILLATIONS 



The equation r/as solved for eight values of The 

roots are given in the following table, the coefficient of 

static stability being evaluated by taking degrees for the 
unit angle. 





V&lues 

Of [X 






V alue s 

Short -period 
oscillation 


of 


A 


Long -peri od 
osci llati on 


n 


* 0.008 




A 


= -3.75 ±3,084 




A 




-0.0408 ±0,2381 




= 0.006 




A 


■ -3.75 ±2.51 




A 




-0.0412 ±0.2131 


n 


= 0.004 




A 


= -3,74 ±1,861 




A 




-0.04 37 ±0.19i 




= 0,002 




A 


» -3.73 ±0.9281 




A 




-0.0508 ±0.1411 


l» 


= 0,001 


A 




-4.58 A = -2.86 




A 




-0.056 ±0.092i 




= 0.0 


A 




-5.18 A = -2.255 


A 


a -0 


.1275 A = 0 


M> 


= -0,001 


A 




-5.6 A = -1.8 


A 


= -0 


.2307 A = +0.075 




= -0.002 


A 




-5.97 A = -1.455 






-0 


.2895 A = +0.127 



The notion studied is stable when the roots are nega- 
tive or v;hen, if imaginary, their real parts are negative. 

The factor e a * or e^* approaches, in this case, zero 
as t increases. 



The time after rhich the factor e at attains the 
value 2 is given by: 

, ln 2 0,692 

i | = = — 

|a| fa I 

TJhen the coefficient of static stability has a suf- 
ficiently large value the characteristic equation defines 
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in general two very different types of oscillation. One 
oscillation is of short period and is very damped; the sec 
ond has a long period and is less damped. As the static 
stability decreases the motion ceases to he oscillatory. 
In the example given, the rapid oscillation loses its char 
actor somewhat earlier than the slow oscillation. The pe~ 



an d da 


mping are 


given in 


the tat>lc 


"be low 




Sho rt - 


period 


long-period 




oscillation 


o s c i 1 1 a 


t i on 




( SO CO 


nds ) 


( seconds ) 


n 


T 


■fi 


T • 




0.008 


2,08 


0.185 


26.4 


17 


.006 


2,51 


.185 


29.5 


16.8 


.004 


3.38 


.185 


33 


15.9 


,002 


6.7 


.186 


44.5 


13.6 


,001 


aperiodic 


68 


12.4 


.00 


ap e r i o d i c 


aperiodic 



TThen [l » 0, the two kinds of oscillations vanish, th 
roots becoming real and one of them zero. 

In passing to the condition of static instability the 
total motion is the sum of four periodic motions, three of 
them corresponding to disturbances decreasing with tine, 
the fourth corresponding to a disturbance increasing with 
time and indicating dynamic instability. The table of 
roots shov:s clearly that the unstable motion corresponds 
to one of the components of the slow oscillation. The 
aperiodic motions which replace the rapid oscillation do 
not cease to be stable within the limits of static insta- 
bility con sidered. 

T7hen the motions are aperiodic the equations defining 
the motion cannot bo put in the sinusoidal form Till a&d 
must renain in the exponential form IT. 
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CHARACTER OF THE OSCILLATIONS 



The character of the two typos of oscillations is 
T7ell knorm. The rapid oscillation is primarily that about 
the center of gravity of the airplane, rzhile the slo^ 00* 
dilation, on the contrary, is one of the path described 
"by the center of gravity in the vertical plane. The latter 
type of oscillation originates in the irregularities pro- 
duced by any disturbance of the equilibrium of the forces 
applied to the airplane* It is possible by approximate 
methods to study each of these types of oscillations inde- 
pendently • 

^l_Sa^id_o^cillat ion.- In studying the rapid oscilla- 
tion: by the method utilized by Munk - taking account of 
the loss in altitude of the airplane - it v;ill be found 
that the restoring moment is not proportional to dC?.j/ di , 
but to 

dO M £C^ dC_^ SM,» gl t 1^ 
di + di di 1 SI V s C z 

and the damping is not proportional to 

dC' z s n ■ 2 
di" 7 Sr^~ 

but to 

dC'_ SU f2 dC 7 
z + — 2 

di 1 Sr 2 di 

This explains why, so long as the maximum lift is not 
attained, the rapid oscillation could be stable even if the 
coefficient of static stability is negative. 

The rapid oscillation, at small angles of attack, is 
strongly damped and the function ? a *p sin (bt + cp) , as a 
result of the values of a and b, respectively, practi- 
cally ceases to have an oscillatory character. If t for 
example : 

a - - 3.75 b = 1.5 6 



then 
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T J a 0^692 s 0#18g 
la | 



which means that after about i/20 of a period tho ampli- 
tude ir, again affected "by a factor 0.5. In the case where 

cp m f the amplitudes start out from p, out decrease 

more rapidly than according to the sinusoidal curve, the 
damped motion "being that shoun in figure 4. The airplane 
is, in fact | very energetically brought hack to the angle 
of attack which corresponds to a zero moment about the 
center of gravity. 

The sane doer, not apply at the large angles of attack. 
At the maximum lift ™- = 0, and tests on numerous models 



show that dC f z /di 1 is also very small at this instant, 
whereas dOjj/di , on the contrary, maintains its former 
value. The theory thus predicts that the rapid oscillation 
may not be damped at the large angles of attack. 

b) Slow oqcHlfrtlqa .** The slow oscillation may he stud- 
ied separately "by examining the modifications undergone by 
the path when the velocity is subject to a disturbance. It 
Will readily be found that the flight path frill be of an 
oscillatory character. These oscillations are necessarily 
accompanied by disturbances of the other variables. Tak- 
ing account of the fact that the airplane strongly responds 
to the applied moments II and of the tendency, through 
the short-period oscillations, to assume rapidly the angle 
of attack of equilibrium, it might be supposed that these 
long-poriod oscillations are effected at a strictly con- 
stant angle of attack. This conclusion would, however, be 
premature. The moment M is, in fact, a function of two 
variables w and q, and we may write: 



The rapid oscillation does not permit the continued 
existence of moments M and therefore, dM a 0. The ir- 
regularity of the flight path, however, corresponds to 
the existence of angular pitching velocities dq , so that 
the condition dM a 0 implies the existence of disturb- 



di 



dM = 



dM 
3w 



dw + 



3M 
dq 
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ancer, of angle of attrck 6w or 5i connected with 6q 

by 

— a 6w + S Sq = 0 

9w 3q 



We thus have : 



I S»t * dC t \ 3G M 

6 q I g- ) = - S w — ~ 

V SI 2 di* / 9w 



3C M 1 
di V 



whence 

30 1.: 

_ _ 1 si 8 IT" 

This condition will only "be true provided that: 

1) p 3 and p 2 o.re in the previously given ratio, and 

2) Sq and 8w corresponding to the short-period os- 

cillation are out of phase by tt. 

In the case of the numerical examples, tire should have: 

p 3 = 0.01385 p 2 for the airplane for which M» = 0,002 
p = 0.0352 p M» = 0» 00 8 

3 2 

It may ho seen from the numerical tables computed, 
that this relation is practically verified and that the 
phase displacements are by tt i 3°, 



LETERMI1TAT 101! OF THE 1:0 T I Oil AFTER INITIAL DISTURBANCE 



In the general motion each of the elementary motions 
is involved to a degree that depends on the nature of the 
initial disturbance. Preliminary studies on automatic 
stabilizers have led us to investigate how a given initial 
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disturbance is distributed between the two types of oscil- 
lations end we here give the results of these computations. 
The ganeral method indicated above has been applied to a 
sufficiently large number of particular cases. We have, 
however, limited our computations, as far as longitudinal 
motion is concerned, to six airplanes instead of eight, 
eliminating the intermediate case p = 0.001 as well as 
that for the unstable airplane \x = - 0.002* The method 
used becomes in fact essentially inapplicable when^the 
disturbances are no longer small. Since certain disturb- 
ances increase when there is instability the method is at 
fault in cases of instability, and we have therefore lim- 
ited ourselves to the study of the unstable case \i - -i 
0.001, v;hich leads to a less rapid increase in the dis- 
turbances . 

T7c shall examine the motion of the airplanes for the 
following initial disturbances: 

1. Increase in the angle of attack caused by a sud- 

den pitching motion. 

2. Increase in the angle of attack caused by a ver- 

tical gust. 

3. Increase in the relative velocity due to a gust 

acting in the direction of the flight path. 

4 # Sudden engine failure. 

In a general way we assume the initial disturbance 
to be sufficiently large. It may be objected that we are 
here departing from the assumption of small motions. The 
size of the disturbance has been chosen, however, with a 
view toward obtaining variables that may readily be handled 
numerically. Since all the equations are linear the ini- 
tial disturbance may be divided by an arbitrary number and 
all the amplitudes will be divided by the same number. 

DISPLACEMENT OF THE AIRPLANE IS UH DISTURBED AIR 



Suppose the airplane has undergone a sudden disturb- 
ance displacing the airplane by a certain angle with re- 
spect to the horizontal and the flight path. We shall as- 
sume a di si)lacement of 0.2 radian in a direction to nose 
the airplane up. The airplane will then undergo a dis- 
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tW*ba&ce in altitude 88 = - 0.2 radian or - 11.46°, and 
a disturbance in the angle of attack U = + 0.2 radian, 
corresponding to 8w = - 0.2 7 - - 8 m/s (18 m.p.h.). 

The numerical values found for p , cp , etc, are given 
in table I, and the diagram giving the disturbances 6u 
of the velocity u, (which nay he assumed to "be the same 
V), 6 6 of the inclination 6, and 8i of the angle 
of attack i, is shown on figure 5. These curves, as 
well as the figures of the numerical table, show how the 
initial disturbance affecting the angle of attack and the 
attitude, is distributed between the two oscillations.* 
The airplane is suddenly raised as a result of the excess 
of lift and tends to nose down if stable - these two phe- 
nomena decreasing the angle of attack. 

The most stable airplane noses down energetically by 
the action of the rapid oscillation. After a fraction of 
a second the disturbance of the angle of attack is prac- 
tically annulled by the joint effect of this diving ac- 
tion and the vertical acceleration which turns the flight 
path upward. 

The strong curvature of the 80 curve corresponds to 
the vanishing of the rapid oscillations. At the instant 
when these vanish the airplane is at an angle of attack 
differing little from the normal angle of attack but 
stalled by 6° and on a rising flight path. A lack of equi- 
librium will be felt in the forces and the long-period 
oscillation will arise from this fact. 

The airplane only slightly stable likewise noses down 
through the effect of the short-period oscillation but the 
action is loss accentuated. The flight path meanwhile 
curves upward. At the instant when the rapid oscillation 
ceases to be felt, the airplane has nosed down only 2°, 
but its flight path has been raised and the airplane will 
be found on a rising flight path of about 9°, the dis- 
turbance in the angle of attack always remaining small. 
The nonequilibrium of the forces is greater than in the 
preceding case, and the resulting long-period oscillation 
will be of greater amplitude. 



On this figure, as well as on those following, the nega- 
tive values of 88 are plotted above the x-axis so as to 
facilitate reading the diagrams. Upwardly inclined mo- 
tions will then correspond to rising curves and vice versa. 
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The long-period oscillation in very clearly brought 
out by the disturbances in the velocity ?, to which it 
corresponds. The disturbances in the angle of attach 
v/hich accompany the long-period oscillation are greater 
on a less stable than on a very stable machine • The case 
of an airplane statically neutral requires no remarks. 
In the case of the statically unstable airplane the dis- 
turbance of angle of attack decreases at first as a result 
of the upward curvature of the path; the latter is not 
stable, however. The airplane enters a condition of con- 
stantly increasing angle of attack, leading inevitably to 
such disturbances that the method of snail notions, after 
a certain tine, ceases to apply. 



DISTURBANCE OF ANGULAR VELOCITY 



The disturbance above considered, suddenly increasing 
the angle of attack and the orientation in space of an im- 
portant angle, is not actually realized in practice. It 
may, of course, be imagined that a localized gust strikes 
the rear of the airplane, but the latter cannot instantly 
attain the final angle of attack assumed. The airplane 
will pass through all intermediate angles of attack and 
the effects of excess lift and of the static stability 
will make themselves felt during these states* 

In order to analyze the caso, let us imagine ^ that the 
disturbance applied is an impulsive angular velocity tend- 
ing to turn the nose up 5q < 0. From the diagrams of V, 
B t i, v:e have found that for 8q = - 0,96 radian per sec- 
ond (r;hich is a considerable value), the ultimate effects 
of the disturbance wore the same as those of the initial 
disturbances of attitude and angle of attack studied above. 
The immediate effects cf the disturbance evidently differ. 
The angle of attack starts from zero and tends to increase 
but the increase is rapidly stopped and the long-period 
phenomena are the only ones that subsist after about two 
seconds. It does not appear necessary to comment any fur-* 
thor on the diagrams of~ figure 6, constructed for the ini- 
tial disturbance 8q ■ - 0.9 6 radian. 

DISTURBANCES OF THE SURROUITDIHG MEDIUM 

The effect of the disturbances in the surrounding me- 
dium may be studied by the following device. Let us as- 
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sume that the ambient air, instead of feeing undi s turned , 
is subject to a velocity U f , having the components u 1 
and w r . Let U of components u and v: he the absolut 
velocity of the airplane. The relative velocity may he 
written : 

and the angle of attack as 

4 ^ W *"» W 1 

u - u T 

It is always possible, in the case of a given air- 
plane, to investigate the mo-tion corresponding to an equi- 
librium of forces in an atmosphere having any absolute 
velocity whatever, since the relative velocity and the an- 
gle of attack are independent of the disturbance velocity. 
It is therefore simple to investigate the motions which 
f airplane would assume during disturbances of the sur- 

rounding medium. 

Lot R be the condition of the airplane for which 
u } x = 0 w ! i a 0. This condition corresponds to the val- 
ues u x and w x of the absolute velocity of the airplane 
In a disturbed atnosphere A 2 for which u r 2 Q n r 2 ^ 0 
the airplane, flying at the same angle of attack and sane 
relative velocity, will be in a state R 2 characterized 
by different absolute velocities u 2 and w 2 determined 
by u 2 s u x * u ! 2 T7 2 a w g * w T 2 . If the airplane flying 

in state R 1 suddenly passes from the atmosphere A x to 
the atnosphere A 2 , it will cease to be in equilibrium. 
It will not be able to maintain its state R x and will 
tend to assume state R 2 . The motions of the airplane 
could be calculated by considering the state R 2 as that 
of equilibrium, and the airplane as deviating from it by 

6u = - u ! 2 
6V7 = - w T s 

This artifice enables us to study the effect of sudden 
gusts on the airplane. 
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ASCENDING GUST 



Let an ascending gust, characterized "by u f 2 = 0 
w r a = 8 m/ s act suddenly on an airplane in horizontal 
flight. The state R x in the atmosphere A x assumed to 
"be undisturbed, is characterized "by a system of given val- 
ues u x , tr Xl 6x (fig. ■?)• The state H x in the atmos- 
phere A 2 is a transitory state and not one of equilib- 
rium. The angle of attack is very large (fig. 8) and the 
flight path of the airplane will he modified. The air- 
plane will try to adjust itself to a new state of equilib- 
rium R 2 in the atmosphere A 2 . It is easy to see that 
this state will be characterized by the same values of the 
relative velocity, angle of attack, and inclination to the 
horizontal as state H x , hut the absolute velocity U a 
will be different. The latter becomes ascending, the air- 
plane being subject to the motion of the surrounding me- 
dium (fig. 9). 

We may study this transitory period by seeking the 
effect of an initial disturbance Sw = - 8 m/ s (fig. 10). 
The airplane will be lifted suddenly and if it is statical- 
ly stable, it will start a diving motion. These two ^mo- 
tions have the effect of tending to close the angle in- 
cluded between the directions OX and V of the figure, each 
of these axes approaching the other. 

If the airplane has great static stability, it will 
oscillate more rapidly and the angle of attack of the 
state will be regained by means of a displacement 8 8 
larger than would be the case if the airplane had been less 
stable. The static stability increases the disturbance in 
the attitude 69 due to the rapid oscillation produced by 
a vertical gust. The axis of the airplane should, howev- 
er, in its final state, regain its initial attitude 6. A 
large static stability therefore results in momentarily 
removing the airplane from its final attitude and increas- 
ing the nonequilibrium of the forces acting on the moving 
center of gravity. This explains why the amplitude of the 
long-period oscillations will he greater on a very stable 
airplane than on a less stable one. A neutral airplane 
will be lifted without its axis OX undergoing any change 
in attitude 86. It is the U x axis which will be effec- 
tive in destroying the disturbance of the angle of attack 
Si. The case of an unstable airplane is sufficiently ex- 
plained by the figure and requires no remarks. 
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HORIZONTAL GUST 



Sy the same procedure we have calculated the effect 
of a horizontal gust of 10 m/s (22 m.p.h.) acting in the 
opposite direction to the motion of the airplane. The 
initial disturbance is an increase 6u. After the effect 
of this disturbance has disappeared the airplane regains 
its horizontal flight at the original angle of attack "but 
at a different altitude. The airplane undergoes at first 
a vertical acceleration. Its angle of attack therefore 
decreases, and if it is stable it noses up and then starts 
the notions characteristic of the long-period oscillation. 
Figure 11 shows the action on the four stable airplanes 
considered, on a neutral airplane, and on an unstable air- 
plane. The disturbance of the horizontal velocity gives 
rise to disturbances in the attitude associated with the 
slow oscillation and their effects in certain cases be- 
come far from negligible. 



DISTURBANCE OF TEE FLIGHT PATH 



We shall consider the case of a disturbance 86 act- 
ing alone - the velocity V , angle of attack i, and an- 
gular velocity q not being disturbed. Such a disturb- 
ance occurs under special circumstances* It would result, 
for example, from the same angular displacement of the air- 
plane and its flight path. The case where an initial state 
and final state are characterized by sensibly the same an- 
gle of attack and same velocity V, but where the paths 
and the orientation of the OX axis differ by the same 
quantity, is one that may occur in practice, namely, that 
of the sudden failure of the engine. Denote the propel- 
ler advance diameter ratio V/nD by Y. On some air- 
planes, when 




the sudden failure of the engine does not change the mo- 
ment about the center of gravity. When this is true, the 
airplane, for the same deflection of the elevator, will be 
in equilibrium with its engine cut off and on a descend- 
ing flight path over which it will travel with the same 
velocity and at the same angle of attack as on the hori- 
zontal path with the engine running. 
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The horizontally flying airplane, subject to sudden 
engine failure, is then at an angle of attack and velocity 
of equilibrium corresponding to the now conditions, out 0:1 
a path not sufficiently descending, and departs from its 
state of equilibrium by an angle which v/e consider as the 
initial disturbance 56 < 0. It is therefore sufficient 
if wo nako the preceding calculation taking 8 6 as the 
di s turban cc . 

Our figures have been drawn for 86 » <* 0.2 radian 
in order to maintain the sane magnitude of initial dia- 
turbance and facilitate comparison with the preceding fig- 
ures. This number does not correspond to any actual value 
of the 'minimum' drag-lift ratio because the latter would bo 
tan 11.4° or 0.225, that is, a very poor drag-lift ratio. 
The amplitudes read off on the diagram should be dimin- 
ished in this ratio if the effect of the engine failure is 
truly to be taken into account. The diagram is neverthe- 
less of some use. It shows that with a nonequilibrium of 
the power as the initial cause the phenomenon is exclusive- 
ly one of slow oscillation. 



AiiaULAH VELOCITIES A2TD ACCELERATION 

The preceding diagrams are concerned with the values 
of the variables V or u, 6, and i or w. tfe have not in- 
dicated the angular velocities, the latter being obtained 
by determining the slopes of the 6 curves. It is equally 
possible to determine the accelerations or the intensity of 
the apparent weight. As regards the sensations of the pas- 
sengers, the two most important variables arc: the incli- 
nation to the horizontal 6 and the magnitude of the ap- 
parent weight J z . 

The quantity J 2 measured by the acccl eromct er , is 
the function 

* < - J z - |f - + e cos e ' 

Under nornr.l conditions J z has the value g and nay 
ccmally well tjc called the acceleration. It can tc calcu- 
lated as a function of the values a, "b , P, and cp given . in 
the taolcs, observing that 
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- P o at [a sin ("bt + cp 2 ) + b cos ("bt + cp g ) ] 
+ p' 2 e a,t [ a « sin (b't + cp' ) + b» oos("b*t+ epO] 



(x) 



It nay also "be determined from the i and V curves 
We have, in fact, 



whence 



and fron the angle of attack i and the velocity, we nay 
obtain F z as 



aV 



2 



whence 



z z 2g 

S C z a V 1 
J = g - 

z ° ? 2^ 



The second procedure for determining J z the more rapi 

if v:c are given the curves of i and V. 

The most important accelerations arc those associated 
with the short-period oscillation. The latter nay he 
found either hy examination of the expression (X) or of 
the curves of velocity and angle of attack. In the ex- 
pression (X) at the start of the phenomenon when t is 
very small, the first term is the pr eponderant one since 
the quantities p, a, and "b entering as factors, are 
larger when they are relative to the rapid oscillation 
than vrhen they characterize the slow oscillation. The dia 
grams also show us that the large variation of i leading 
to large values of C z and J z arc due to the rapid os- 
cillations. In the long-period oscillation the angles of 
attack and velocities V vary simultaneously. Since the 
velocity V increases when the angle of attack i de- 
creases, the variations in the acceleration arc not very 
large. It is important to note that the acceleration de- 
tcrnincd by the preceding ncthod, in the case of a verti- 
cal gust, corresponds to the roughest sort of calculation. 
It is assumed that the gust is suddenly set up, that the 
airplane iron the first instant lies entirely within this 
gust, and that any phenomenon of elasticity of the wings 



IT.A.C.A. Technical Llcnorandun Ho. 867 



23 



Trill not dininish or complicate the action. Numerous 
works "based on more probable assunptionc have "been pub- 
lished. In a remarkable paper published many years ago, 
S. B. TTilson assumed as the law for the set-up of the 
gust, the following expression: 

w = J (1 - e^ rt ) 

Other authors have assumed a sunusoidal law which is gen- 
erally arbitrary. The latter assumption leads immediately 
to resonance phenomena if the period of the airplane is 
the same as that of the gust. Fisher, Bryan, and Jones 
have studied the effect of sudden and gradual gusts by 
methods which have become more and more accurate. Kussner, 
in a very well-known work of his, takes into effect the 
elasticity of the wings. Since our object is to compare 
the effect of different disturbances, we have not gone 
into the subject in more detail. 

EFFECT OF THE ACTION OF THE CONTROLS 



The method outlined above enables us in certain par- 
ticular cases to determine how an airplane responds to a 
manipulation of the controls. Let us examine figure 13, 
which gives 0 M as a function of i and p. Assume the 
airplane to be flying level at the angle of attack ± 1 
with the control deflected at angle 3 X , this condition 
corresponding to a velocity 7 1 . If the pilot gives the 
control a deflection 3 2 , the point A passes to B. At 
B the moment M is no longer zero and the point will be 
displaced from B to C along the curve 0 2 . The air- 
plane will be in rotational equilibrium about the center 
of gravity only at the angle of attack corresponding to^ 
point C. Let i 2 be this angle of attack, to which will 
correspond a speed V 2 different from T 1 . We shall as- 
sume that the maneuver is made without altering the throt- 
tle and that the velocity l l is greater than the veloc- 
ity corresponding to minimum power. Under these condi- 
tions the useful power sufficient for making the airplane 
fly at the velocity V x will not be sufficient to make 
the airplane fly at velocity 7 a . The path of the airplane 

will not be able to remain horizontal but will incline by 
the angle £. The axis of the airplane will be lowered by 
£ + (± 1 - i 2 ). At the instant when the pilot in deflect- 



24 



IT.A.C.A. Technical Memorandum No. 867 



ing the control to 0 a , passos fron point A to point B 
he may do considered to have departed from his final state 
C by amounts: 5i = i 1 -i 2 , 8T m Y x — 7 a , 56 = (i x - 
l a ) + ti wh ere £ is computed from the polar of the air- 
plane and the variation of the thrust T with velocity. 

We have applied this computation to a numerical exam- 
ple, choosing the same airplane as above. At about the 
velocity of 40 m/ s (90 n.p.h.), it was found that states 
differing from each other by a difference in velocity 5V = 
±1 n/ s are characterized by a difference in angle of at- 
tack of ?0, 308° . 

Similarly, basing our computations on an assumed po- 
lar, wo found that a variation in velocity of ±1 m/s re- 
quires that the slope of the flight path vary by +0.0070, 
corresponding to d| '« - 0.400°. These data were applied 
by us to the numerical example by multiplying them by the 
factor necessary to obtain amplitudes that may easily be 
studi ed . 

Consider, for example, the pas sago from the level 
flight condition at velocity V = 27 m/s to a state where 
V = 40 n/ s by means of a deflection AS required to pro- 
duce a Ai = - 4°. If this deflection A3 is suddenly 
applied, the airplane is deviated from its final condition 
by: 

67 « - !«; fii = + 4°; || = - 5.2° 

Nov/ 8© d $1 - gi f s0 that 66 « - 9.2°, one portion of 
which represents the change in angle of attack, and the 
other the change in the flight path. It suffices to deduce 
the amplitudes of the disturbances of figures 10, 11, and 
12 in tho desired ratio and add them. Figure 14 (continu- 
ous line) shows the result of this computation for the air- 
plane of coefficient \x = 0.004. The diagram of the angles 
of attack shows that the final angle of attack is not es- 
tablished rapidly. TCo may trace separately the effect of 
each of the elements of the ini t JLal di s turbance on the vari- 
able ^5i (fig. 15). The effect of the initial disturb- 
ance in the angle of attack decreases very rapidly. Tho 
initial disturbances of velocity and attitude* however, also 
have an effect- on the variable i - an insufficient veloc- 
ity tending to increase the angle of attack, and similarly 
for a not sufficiently descending path. These two initial 
disturbances act on the long-period oscillation. It may be 
seen, in the example chosen, that in spite of the rapidity 
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with which the short-period oscillation tends to bring the 
airplane to its final angle of attack, the airplane may in 
certain cases remain considerably deviated from its final 
position by the effect of the components of Si that depend 
on the long-period oscillation. 

Remarks 

1. If we assume that the deflection Ap has not been 
instantaneously applied, but in a progressive manner (for 
example, one-third" at time zero, the second third at time 
0.5 second, and the remainder at time t ^ 1 second), the 
curve of variation of angle of attack loses the abrupt char- 
acter which it has during the first second and becomes reg- 
ular with. its appearance not appreciably modified there- 
after. This condition corresponds to the dash-dot curve of 
figure 14. 

2. The airplane will more quickly attain its final 
angle of attack if the pilot, to start the maneuver, gives 
the control a greater deflection than the amount A0, which 
he must maintain at the end of the phenomenon. 

3. The computation is made on the assumption that the 
pilot does not vary the engine throttle. If the pilot in- 
creases the power the final path will no longer be descend- 
ing* The effect of such a maneuver would have to be sepa- 
rately studied. 

1.10 RE COMPLICATED MAKIPULATI01TS OF TEE CONTROLS 



The same procedure may, of course, be applied to ana- 
lyze the effect of more complicated maneuvers. Let an 
initial state be characterised by a given velocity attitude, 
angle of attack, and deflection (3. For example, in tne 
case of the airplane considered above, V = 27 m/ s ; i = + 
4°; 8 -a - 4°. 

The pilot deflects his control by A3, and the air- 
plane starts the motions which tend to bring it to a final 
state B. He applies this deflection p + A0, however, 
only within a limited time interval, and re-establishes the 
deflection 3 after n seconds. At this instant the^ air- 
plane has not yet attained its state 3 but is in an in- 
termediate state. We see on the figure that the momentary^ 
characteristics after 7 seconds arc: V a 34 #7 m/sj i = 1 ; 
6 es + 13.4°; q = + 1° per second. 
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The final state of equilibrium, however, defined by 
the position of the control, has "become state A, and the 
airplane deviates from it "by 6T - + 7,7; 8i = - 3°; 
86 = 17,4°; Sq = + 1° per second. Taking these values as 
an initial disturbance with respect to state A, the 
changes in the variables may "be easily determined. The re- 
sults of the computation are indicated as discontinuous 
lines on figure 14. The computations have been made -neg- 
lecting the effect of the disturbance 8q, which is less 
important than the three preceding ones. 

Remarks 

Hm The curve of the angle of attack drawn as a thin 
dot-dash line corresponds to the assumption of a more grad- 
ual maneuvering of the elevator. 

2» We have calculated the curve of J z , the apparent 
weight or acceleration, for the maneuver described - con- 
sisting of nosing down the airplane for 7 seconds, then 
nosing up. The result of the computation is given on fig- 
ure 16. 

Tests 

It was possible for us to obtain records of the varia- 
bles defining the longitudinal motion of an airplane. The 
test was carried out on a Fairey "Fox" with the Bouny 
equipment, which will be described in detail in a succeed- 
ing Bulletin. Records of "phugoid" oscillations are suf- . 
ficiently numerous. Nevertheless, we believe it useful to 
publish here the results of the measurements we have made. 

Although this test has been conducted on an airplane 
differing from the one employed in our calculations, the 
curves obtained show a qualitative correspondence with the 
calculated curves that appears to justify comparison. The 
polar and the curve of moments were known from a tunnel 
test conducted under the usual conditions on a model not 
provided with a propeller. The airplane itself had a load- 
ing of 58.5 kg/m 2 . At the speed of 45.2 m/s, it flew at a 
lift coefficient of 0.458, corresponding to an angle of at- 
tack of 3.6°. The coefficient of mean static stability 
measured in the interval i = + 6 to i = - 4 may be de- 
termined from figure 17, 
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The coefficient of effectiveness of the elevator for nega- 
tive deflections was: 

-~ = 0.01308 
dp 

The conditions of equilibrium , according to the tunnel q 
tests, bKott that the flight at G 2 = 0,458 and i m 3.6 , 
should he made at an angle of elevator deflection of 3 = 

- 2.55°. Actually, equilibrium was obtained at p = 8»6 
which is a sufficiently good agreement. 

The test consisted in carrying out the following ma- 
neuver: Having attained the steady state at velocity V m 
45.?. m/s, which here is state A, the pilot pushes on the ^ 
stick and lowers the elevator, the mean deflection "be coming 
p s m 1.5°, the pilot not touching the throttle. The air- 
plane dives and tends toward a new position of equilibrium. 
The pilot, however, does not wait until the final state B 
is attained. After having maintained the deflection^ 0 = 

- 1.5° for 7 seconds, he pulls hack on the stick, fixes it 
in its initial position and allows the airplane, after &©me 
oscillations, to return to its initial state A. The test 
is carried out while recording: 

V is the velocity on the flight path. 

V , vertical component of the velocity measured 
v 

with a variometer. 
P f setting of the elevator, 
i, angle of attack, 

J z , component of the apparent weight in the direc- 
tion of the OZ axis. 

J x , component of the apparent weight in the direc- 
tion of the OX axis. 

q, angular velocity about the OY axis. 
We thus have 7 test curves. 

The variometer is an apparatus whose readings show con 
siderahle lag, and can only he used when corrected. The 
curve of corrected vertical velocities is shown on one of 
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the iiagraag of figure 18, A kaowle&ge of Y T and ¥ 
gifes tlie angle of elope of the path and the angle 6, may 
he calculated iron g and i. Finally, having deter- 
mined G, it ir, possible to calculate q and check the 
readings of the indicator of the angular velocity of pitch . 

TCe shall nov; consider the curves distinguishing three 
portions : 

1« 'The flight is assumed rectilinear and uniform at 
state A, the deflection being -»$«6 9 the part 
of the curves to the left of point 1. 

2 # A ivcriod of 7 seconds duration, during which the 
airplane undergoes maneuvering rzhich would lead 
to state 3 if the deflection 0 a? ^ 1.5° -vere 
maintained. This period corresponds to the por- 
tion 1-2 of the curves. 

3. She period following the return to the setting- 

p = - 3 # 6°, during which the airplane oscillate 
end tends to regain its initial state A - that 
is, the part of the curves to the right of point 
2. The relation: 

£a p - G Si! 

g S 2 rg 



permits us to determine C z and hence, i. 

For the angle- of -attack curve the computation gave a 
curve which differed slightly from the one recorded, the 
latter lagging behind the computed curve and at the mini- 
mum angle of a,ttr,ck not soiling down so low. It was possi- 
ble for us to show that there was a systematic error in 
the reading of the vane due to a play of a,bout 1*5 f a 
fact which explains in part the lag of the record. The 
Bouny acceleroncter is an instrument with which we .arc very 
familiar and whose readings arc very accurate. Under these 
conditions we think that the computed curve corresponds to 
that of actual angles of attack. 

The curve cf angular velocities has "been obtained with 
the aid of an apparatus designed to measure the angular ve- 
locities which are produced in spinning and which are of 
the order of 1 to 2 radians oor second. This apparatus 
did not have the desired sensitivity for measuring angular 
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velocities which are of the order of 0*05 radian per second 
We were able, however, to calculate the angular velocities 
q from the angles G. The angular velocities thus calcu- 
lated are, however, obtained after a rather larger number 
of operations, among which is the correction for the indi- 
cations of the variometer. Under these conditions, the 
agreement of the two curves is not to he considered as very- 
poor. 

We shall now consider the curves in detail. After the 
elevator has "been deflected from P ■ - 3#2° to (3 = ~ 1*5 
the airplane tends toward a final state B, The latter, cor 
responding to a displacement A3 a 1*7° , is characterized 
"by a final angle of attack: 

i f rz 3.6° + Ai 

where 

Ai = - .Otf&gg A3 = - 5° 
0. 00477 

Therefore, l f = 3. 6° - 5° = 1.4° 

corresponding to a lift coefficient: 



C z = 0.458 - 0.0695 Ai = 0.111 

and a velocity V a 92 m/ s 

The airplane used is thus extremely sensitive to the 
elevator controls. The deflection was applied during 7 
seconds, within which time the airplane "began a series of 
motions which would have led to state 3 if the deflection 
had be e& applied long enough. This period is characterised 
"by the immediate appearance of a component of the velocity 
directed downward "by the decrease in the angles of attack 
and acceleration J* z , and by the increose in the velocity 
along the flight path. We nay note that the increase in 
the velocity is made at first rather slowly, thus confirm- 
ing the computations given in figure 14. 

Let us now examine the effect of the return to the in- 
itial elevator deflection. At the instant when the stick 
is thrown hack to its original position (point 2), the 
curve of accelerations is suddenly modified, this bein£ the 
curve which indicates most exactly the instant when the 
maneuver has "been executed. The angle of attack is instant 
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ly changed, the velocity V continues to increase for sev- 
eral seconds, and the airplane continues with a vertical 
velocity of descent - the latter, however, attaining its 
maximum a little after the maneuver. The different varia- 
bles vary in such a manner as to regain the values corre- 
sponding to the initial state. The diagram shows up the 
long-period oscillations with particular clearness, a large 
amplitude "being obtained, due to the size of the initial 
disturbance imparted to the airplane. 

T7c do net consider it necessary to comment further on 
the experimental curve obtained. The curve has the sane 
appearance as that computed and drawn on figure 14 for a 
similar maneuver but corresponding to a different airplane, 
and thus proves the value of the computations for the theo- 
retical determination of the flight paths. 

It is part of our program to carry out complete compu- 
tations, using the characteristics of the airplanes avail- 
able for our tests. 



III. LATERAL OR TRANSVERSE MOTION 



Method and Equations 



The equations of translation along the OY axis, and 
rotation about the OX and OZ axes are those that determine 
the lateral stability. These equations are: 



(ft + rtL ~ pw ) = Y ~ P Gin ^ 



g Vdt 



A || - E — + q r (C - B ) - Epq = SL 



C _ ~ E S£ + pq (b - A) + Eqr = EH 
dt dt 

where A, B, and C are the moments of inertia. If the 
axes of coordinates coincide with the principal axes of in- 
ertia the product of inertia E is zero, and this we shall 
assume in what follows* 



To the above equations must be added two geometric re- 
lations connecting the angular velocities p and r with 
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the derivatives and and resulting from the dof- 

dt dt 

initios of the rotations: 

p = if CC3 6 - sin 6 cog « 
dt dt 

r = 4? sin B + cos 6 cos qp 
dt dt 

These relations may he written: 



dc? 

--'r = P cos 6 + r sin 6 
dt 

M = ( r cos 6 ~ ? sin 6 ) 

dt cos cp 

The system of equations is of the form: 

f x (v, p, r, co, 



dv 
dt 



dn> 
dt 

dr 
dt 

dp 
At 



f 2 (v, p, r, cp, 40 
f 3 (v, py r, q> f • 
f 4 ( v , p , r \ <p \ y ) 



^ - (t, Pi CD, \|/) 



dt 



"5 



and may he linearized as in the study of the longitudinal 
stability, the disturbances 6v, Bp, Br, Sep, 8\$f be coming 
t he var iahl o s . 

Etie integrated system, after the above equations have 
been reduced to a linear system, depends on an algebraic 
equation of the fifth degree in A in place of one of the 
fourth degree. This equation will ho: 
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dfx 




sf i_ 

*L- 




6v 


ap 




Sep 


3\|/ 


3f a 


— - - A 


af P 


af P 

a 


_.a. 


3v 


ap 


ar 


9cp 


3\|/ 






af 3 

--5- _ A 


3*3 


a*. 


av 


ap 


3r 




a\|/ 






3f 4 


Hi - A 


3f 4 


av 


ap~ 


3r~ 


Bqp 


d\|/ 


3^5 


af 5 


9*8 






av 


ap 


3r 


3cp 





= 0 



It will "be found, ho\:evor , that the derivatives of the 
functions with respect to $ are sero. The fifth degree 
equation admits of a zero root X n G« This facilitates 
the analytic study since the characteristic equation "be- 
comes one of the fourth degree when this particular solu- 
tion is eliminated f and the mathematical study is made by 
methods similar to those employed in the study of the lon- 
gitudinal stability • 



The existence of this particular root corresponds to a 
certain physical fact. In the study of the longitudinal 
motion it is found that certain projections of the external 
forces depend on the angle 6. When 0 is not zero, the 
axis of the airplane is inclined upward or downward with 
respect to the horizon, and the weight has along the OX 
axis a component which adds to or subtracts from the pro- 
peller thrust. When an airplane is dynamically static, it 
returns to its initial state after a series of oscillations. 
The forces acting on the airplane should therefore regain 
their initial values, and this result can he obtained only 
if the airplane regains its initial longitudinal attitude. 



In the study of the lateral motion we meet with the 
angular magnitudes cp and v|/ . The component of the forces 
along the OY axis depends on the angle cp, since the weight 
has a lateral component when the airplane is inclined. If 
the airplane is dynamically stable it returns, after any 
disturbance, to its initial state and should therefore re- 
gain its initial inclination cp . The same docs not hold 
for the angle \|/ . Whatever the final position of the air- 
plane the projections of the weight on the axes are in&e- 
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pendent of the rotation njr. There exists uo force or mo- 
ment which is a function of the azimuth. Kcncc, when an 
airplane | dynamically stable, returns after a disturbance 
to its initial state, it should regain a notion character- 
ized by the same velocities and by the same angles p and 
6 as the initial motion out not necessarily by the same 
angle \|/ , The airplane, after a disturhar.ee, does not pos- 
sess any weathercock stability ~ all the derivatives with 
respect to $ of the moments and forces "being sere. The 
existence of a solution A = 0 is the mathematical conse- 
quence of this fact. 

The Integral system will he analogous to that obtained 
for the longitudinal stability and will bo written: 



= 




Ci 


gXjt 






c a 


c ^s" 0 






eM 


+ 




r\ 


c _ ^4t 


Sp = 


K 


Ci 


iM 










* N 


c 3 


G *3* 


+ 




c 4 


e ^4* 


69 m 




Ci 


e 








c >. 2 t 


+ n 3 




C 


J- 
1 






e ^4* 


Sep m 




r\ 


e ^ 




r -2 




S 3 


+ H 


s 3 


e A 3 t 




n 4 


c 4 


e x 4 t 



The disturbance 6\|> can be obtained only by integration: 

It is known that the characteristic equation in A al« 
ways admits in the case of lateral stability a pair of im- 
aginary roots A x j2 and two real roots ^3 <and 

motion will result from the su.perpo si t i on of an oscilla- 
tion and two aperiodic motions. 

STATIC STABILITY IH LATERAL EQUILIBRIUM 

The static stability is a phenomenon by whidi a dy- 
namic restoring moment arises when the airplane has beer- 
sub ject to any displacement with respect to its flight 
path. The displacement as regards the lateral motion shows 
up by an angle of yaw j. The ya\7 takes place whenever 
the aerodynamic velocity V departs from the plane of sym- 
metry of the airplane and gives rise to two moments L and 
K . 17c shall take: 
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L = C L Sb g- 

V 2 

N = C ¥ St ^~ 

where the reference length b is by convention the span 
of the airplane. 

The derivatives ~- and fS determine the static 

oj dj 

stability of the airplane. They are the coefficients of 
static stability of roll and of yaw . Since j = v/V, we 
have al so : 

aG|, b i a Cl 
a^ V sj 

ac 

It is readily seen that is positive when the cf- 

dj 

feet of the lateral surfaces of the airplane situated above 
the OX axis in greater than that of the lateral surfaces 

below the axis and that — - is positive when the effect 

of the vertical surfaces situated at the rear is predomi- 
nant. The forn of the wings enters as an important factor 
in the lateral static stability. The dihedral increases 

a o Tj a o i j 

aj~ T7llilc thc Plan forn affects These coefficients 

are capable of numerical evaluation. Tests show that thc 
moments L and II depend also on the angle of attack, so 
that it must always he specified, when a numerical value 
is given, to which angle of attack it corresponds. 



AIRPLANES STUDIED 



Two sets of airplanes were studied - differing in 
their coefficients of static stability of roll and yaw ~ 
the same airplanes being used as in the study of longitudi- 
nal stability. It is necessary to give the span Tr and it 
will be assumed as equal to 10 meters, while the two radii 
of gyration will be: 
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r 3 s c2 . o m 



The clerivat ivo s of the forces and aerodynamic moments 
are given "by: 

Y. r . tts sb §2 L- r . », g - »; t| 

In our conputations we shall assurao the constant values 

a x = - 0,60 a 5 - - 0%M*S 

a 2 = + 0.013 a 6 = - 0.056 3 

ft ts + 0.065 a 8 0 + 0,0168 

a 9 a - 0.044 

It may readily be seen that : 

9 0 L 3 Or 

a4 „ and a? = _ 

where the angles aro expressed in radians. 

The stability of roll will bp characterised "by the fol- 
lowing values of a : 

a 4 = 0, 0.02, 0.04, 0.05, and 0.08 

In investigating the effect of these variations in stabili- 
ty, we shall comhine it with the constant value a ? = 0.04. 

The directional stability will oe characterized fey the 
following values of a ? : 

a = - 0.02 -0.01 0 +0.02 +0.04 +0.0 65 

7 

In varying a , W6 shall suppose a 4 constant and equal to 
0.04. 
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Re mar h 

TTg assume, in fact, that the coefficients of stability 
may vary without having the quantities a 5 , a 6 , a Q , and 
a g , which define the aerodynamic damping, vary , It is im- 
portant to note this assiimption. 



SOLUTION OF THE CHARACTERISTIC EQUATION 



The movements executed "by the ten airplanes studied 
arc not all stable. The criterion of stability of Routh 
requires, when the characteristic equation is put in the 
form (VII), that each of the quantities 

A i, A*, A 3f A 4 , and S = A, r - ^ A 

he greater than zero. 

TCe shall take the coefficients of static stability as 
variables 

a 4 - x a? =y 

The expressions Aj , A 2 , A 3 , A 4 , and R become func- 
tions of x and y. Let us consider the region of the 
plane determined by the extreme values of x and y from 
-0.02 to + 0.10 (fig. 19). 

The quantities A; and A 2 are always > 0 for the 
values of the variables x and y considered. The lines 
Ax a 0 and A 2 = 0 do not cross this region of the plane. 
The other three relations: 

A 3 = 0 A 4 = 0 R = 0 

on the contrary, define lines which cross the region of the 
plane considered. 1c shall determine a region where sta- 
bility exists. 

Plotting the points . characterizing the ten airplanes 
studied, it may be seen that six are dynamically unstable, 
and only four are stable. The roots of the equations and 
the characteristics of the motions are, given in the followin 
t able s : 
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VARIABLE STABILITY OF ROLL 
Value of a 4 Imaginary roots Real roots 



0.08 


\ 

-0. 325 2 


±1.442i 


\ 

-0.02177 


X 4 
-3.97 


.06 


-.35 95 


±1.40 Si 


-.007 


-3.9 3 


.04 


- . 385 


±1. 3661 


+.01045 


-3. 39 


,02 


-.4126 


±1. 310 i 


+.3035 


-3. 35 


.00 


-.441 


±1.2721 


+ .492 


-3. 815 


The period T 


of the o 


sci Hat ion , 


and the time 


V* 



required for reducing each elementary motion "by half its 
value, or doubling it, according to whether the movement! 
are stable or unstable, are given in the table below: 



Value of a. Point Oscillation Aperiodic motion 

T T it T 4 T 4 



a 



2 



0.03 A 4.31 2.13 31.8 0.175 

.06 3 4.45 1.93 99 .1765 

,04 C 4,53 1.8 66.3 « .178 

.02 D 4.78 1.68 2.28 $ .181 

CO 

.00 S 4.9r? 1,57 XM * .181 

V/hen the directional stability is variable the roots 
are as follows J 
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VARIABLE DIRECTI01TAL STABILITY 



Value of a T 


Imaginary roots 


Real 


roots 




x s 


2 






+0.065 


-0. 3903 


±1. 68Si 


+0.0259 


-3.895 


+ . 040 


-,385 


xl. 366i 


+ .01045 


-3, 89 


+ .020 


-.371 


±1.022i 


-.0209 


-3.89 


0 . 


-.264 


±.556i 


-.2317 


-3,89 


-.01 


+ .04265 


± . 326 3. 


-.846 


-3,89 


-.02 


+ .199 


±.2531 


-1.58 


-3.89 



and the movements are characterized by- 
Point Oscillation 



Valtie of a 7 



+0.065 

+ .04 
+ .02 

0 



t i 
" 3 



3.73 1.78 



2 
3 
4 



4,6 
6,15 
11.3 



1,80 
1,87 
2.6 2 



Aperiodic motion 



26.8 2 

CO 

a 

66.3 * 



32,2 
2.91 



T 1 



0.178 

.178 
,178 
.178 



-.01 



-.02 



19.25 16.25 



24 



to 

48 S 



795 



,437 



,178 



,178 



In what follows we shall denote as the oscillatory mo — 
tion that which is determined "by the imaginary roots; spiral 
motion as the aperiodic motion which is defined "by the root 
X ; and damped motion as the aperiodic motion defined "by 

the root X 4 . We find that the notion corresponding to the 
fourth root is always stable and strongly damped, and that 
its characteristics are independent of the coefficients of 
static stability of the airplane. 
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YThcn the airplanes A, 3, 0, L , S arc successively 
examined, it is seen that the oscillatory motion is stable 
in all cases with practically the oa:i-e period and effec- 
tively damped* The aperiodic motion determined by X 3 is 
stable at the beginning "but not to a great extent, since 
the disturbances do not decrease "by half after a sufficient- 
ly long tine. This motion "becomes unstable as soon as the 
limit A 4 - 0 is crossed, the disturbances doubling in 
times which become smaller and sr.ialloi-. This instability 
is the well-known spiral instability. 

Let us consider the airplanes 1, 2, 6 - differing in 
their static directional stability. The aperiodic motion 
X 3 is unstable at the beginning, the characteristic points 
lying outside the region of stability. The airplane is sub- 
ject to spiral instability but the latter disappears as 
soon as the boundary A4 - 0 is crossed and the motion A 3 
becomes very stable. As the directional stability decreases 
the period increases, the damping decreases, and the oscil- 
lation becomes unsta/ole at the instant the line R m 0 is 
crossed. 

CHARACTERISTICS OF THE HOT I OH 



Dynamic considerations thus indicate that the stabili- 
ty coefficients — ~ . cannot be given arbitrary val~ 

uos, A siiAple qualitative explanation of these facts may 

dC L dC^- 

readily be given. If the two conditions > 0 jj-j > 0 

are nerely satisfied without any further investigation of 
the movements which will be set up, it is possible to have 
the airplane dynamically un stable. 

Let us o::aninc the case of an airplane banked to the 
left 8cp > 0. It will then be under the action of the 
weight on that side and a disturbance 5,1 > 0 will take 
place (fig. 20). The latter disturbance will have two ef- 
fect s • 



1. The airplane will have a tendency to right itself 
under the action of -™ , the yaw toward the left tending 
to incline the a,irplane toward the right. 
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2. The airplane will have a tendency to yaw to the 

left since it "behaves like a weather vane, -~ "being Pos- 
itive, oJ 

There thus arises a new disturbance, a yaw 8r > 0. 
In this motion the right wing will he outside of the turn 
and will he displaced more rapidly than the left wing. Its 
lift will he greater and it will tend to raise itself fur- 
ther, thus increasing the lateral hank. The moment L is 

a function of r, the derivative — ~ is negative, and 

the yaw Sr has a tendency to right the outside wing, which 
here is the right wing. Iwo opposite effects are produced 
and, depending on the proportions of the airplane, one or 
the other will predominate. If the effecrt of the static 
stability of roll is greater, the airplane after yawing 
will regain its initial state - flying, however, in a direc- 
tion different from that preceding the initial disturbance 
6cp . If the effect of the angular velocity r on the roll- 
ing moment is greater than the effecrt of the sideslip j , 
the second effect will predominate. The hank of the air- 
plane will increase, the airplane will describe a spiral 
path, and will then be dynamically unstable. In spite of 
its static stability of direction and roll, it possesses 
spiral instability since the first stability is too great 
with respect to the second. 

Similarly, we can explain how the oscillation assumes 
more and more unfavorable characteristics when the stabil- 
ity of roll is too great compared with the yaw stability. 
Let us imagine that the airplane yaws to the left, the OX 
axis of the airplane being oriented toward the right with 
respect to the flight path. If dC L /aj is high, the air- 
plane will incline energetically to the right and the roll 
will be positive. The tendency to yaw to the left will, on 
the contrary, be small since d%/$J is small. This yaw 
will not develop, as a secondary effect, any important nega- 
tive rolling moment. Since nothing opposes the motion of 
positive roll toward the right, the airplane will bank to 
this side but it should then begin to yaw in this direction 
and the same phenomena will be reproduced in the opposite 
sense* The airplane will have a yawing motion on which is 
superposed a continuous lateral swinging. For too small 
values of dC];j/dj, the successive amplitudes increase and 
the oscillation is unstable. 
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shall here discontinue our qualitative explanations 
and return to the results of our computations. The first 
point is to determine the relative importance of the ele- 
mentary notions "by calculating their amplitudes for the 
given initial di sturbancc s . 



G-RAPHS 0? THE MOTIONS 



We shall consider three initial disturbances, namely: 

1. A disturbance 8v of 10 meters per second, corre- 

sponding to an angle of yaw of 0.25 radian. 

2. A disturbance Sep of 0,2 radian. 

3. A disturbance or of 0.2 radian per second. 

It did not appear necessary to calculate the effect ^of 
a disturbance 8p~/ Such a disturbance would cause the air- 
plane to' bank and would produce in a short time a disturb- 
ance 6cp. TCe prefer to take Sep directly as the initial 
motion to be considered rather than study the sequence of 
successive values of 8cp corresponding to an initial 5p. 
The amplitudes of the components of the motion arc given, 
for each initial disturbance, by the tables which follow. 
The resulting motions are shown plotted on figures 21 to 
23. Several component motions are superposed on these 
plots. Examination of the curves leads to the following 
conclusions. In cases where there is stability, for air- 
planes A, B, 3, and 4, the curves do not present any essen- 
tial differences* Spiral instability exists for the air- 
planes 1, 20, E, and D. It is not in evidence for airplane 
2C, since it is very small and the curves are not prolonged 
a sufficient length of time for the instability to show Up* 
The spiral instability is clearly evident on airplanes E 
and D~ Similarly, the oscillatory instability of airplanes 
5 and 6, is clearly indicated. From a direct study of 
these figures, however, the phenomena cannot be analyzed. 
To see more clearly, we shall investigate separately: 

A. How the initial disturbances considered are dis- 

tributed among the different partial component 
notions . 

B, « 'Jhi each partial motion which are the variables most 

affect ed. 
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Tie shall thus determine the relations existing, in each 
component notion, "between the values of the different vari- 
ables* We know that, for each airplane, these relations 
are constant and independent of the initial disturbance. 

A. DISTRIBUTION OF THE INITIAL DISTURBANCES 



Let us examine how the initial disturbances considered 
are distributed among the partial motions. This distribu- 
tion will depend on the initial disturbance considered. 
With the aid of numerical tables, we may easily show how a 
given disturbance Sep is distributed among the amplitudes 
p 4 , n 3 C 3 , and n 4 C 4 . Similarly, we shall investigate how ' an 

initial disturbance 8r is distributed among the ampli- 
tudes p 3 , m 3 C 3 , and m 4 C 4 , and hor; an initial disturbance 

6v is distributed among p f 0 , and 0 4 . These distribu- 
tions will occur in proportions which depend on the aero- 
dynamic characteristics of the airplane considered. 

T7e find that in each case the amplitude of the motion 
corresponding to X 4 is very small. It is represented by 

the hatched region on the figures. This explains the fact 
that the damped motion docs not play any important part 
in the phenomena. The amplitude of the spiral motion cor- 
responding to A 3 is very la.rge when the initial disturb- 
ance is an inclination 8cp. When the initial disturbance, 
on the contrary, is a sideslip v, the oscillatory motion 
is the greater. In the case of an initial disturbance 8r 
there is some spiral motion, but the oscillatory motion 
pr edominat e s . 

B. VARIABLES AFFECTED III EACH KOTIOIT 



For a given airplane, we can characterize each motion 
by the ratios of the values of the different variables at 
the same instant. These ratios are independent: of the in- 
itial conditions and depend only on the stability parame- 
ters of the a,irplano. 

The amplitudes of the disturbances of the variables 
v» P» r » cp associated with the damped motion, are propor- 
tional to the quantities C 4 l 4 C 4 m 4 0 4 and n 4 C, . Those as- 
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80 ci at ed with the spiral notion are proportional to 
c 3 l 3 C 3 :2 3 C 3 ^3 c 3 • Le * us choose as the term of comparison 
the value of the disturbance of lateral inclination nO 
(of desired subscript). The disturbance of lateral veloc- 
ity C gives an angle of yaw when divided by the velocity 
7. The quantities: 

1 JL - 

V nC Vn 

give the constant ratio which connects, in the elementary 
motion considered, the angle of yaw with the angle of lat- 
eral inclination. The quantities l/n give the constant 
ratio between the angular velocity of roll and the angle of 
roll, m/l the angular velocity of yaw to the angular ve- 
locity of roli, m/n the angular velocity of yaw to the 
angle of roll. From the numerical tables we nay derive the 
f oil owing re suits. 

1. Damped 1'otion 

The ratios ™- , do not vary much when the 

7n 4 n 4 n 4 

different airplanes studied are examined. It is not nec- 
essary to follow the changes of these ratios with the co- 
efficients of static stability, and it will be sufficient 
to consider mean values for characterising the phenomena. 
We find: 

- + 0.147 

Vn 4 

In the damped motion a positive bank (to the. right) 
is associated with a positive yaw (to the loft), the angle 
of yaw being always 0.147 times the angle of bank. The 
relation 




shows that a positive bank corresponds to a negative veloc- 
ity of roll tending to make this disturbance disappear 
very rapidly (in 1/3.8 seconds if the angular velocity were 
constant instead of obeying an exponential law). have, 
moreover : 

— * 0.147 and f& » - 0.0384 
n 4 M 
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These ratios show that a negative rolling motion is associ- 
ated with a positive hut very small yawing motion. 

2. Spiral Motion 

For the spiral motion the ratios vary from one airplane 
to another. Let us draw up the tables "below: 



Airplane 


A 


B 


C 


D 


E 


1 


-0 .0536 


-0.0576 


-0.0663 


-1.985 




-1. 


41 


la. 


-0 .0217 


-0.007 


+0.0101 


+ 0.3035 


+0.49 2 


5l 


-0.216 


-0.216 


-0.2315 


-4.17 




-2. 


3 


m 

-3. 

h 


+ 10 


f 30 . 9 


-22.9 


-13.75 




-4. 


6 


Airplane 


1 


2 


r - ■ 

U L__ 


4 


5 


6 


1_ 
Yn 3 


-0.0416 


-0.66 3 


-0 .0115 


-0.42 


-1 .195 


-1.66 


k 

n 3 


+0.0259 


+0.0101 


-0.0209 


-0.23 


-0.846 


-1.58 


2k 

n 3 


-0.222 


-0.2 315 


-0.229 


-0.207 


-0.81 


-0.91 




-8.6 


-22.9 


+10.95 


+0.893 

I— 


+ 0.96 


+ 0.57 



We sec that the ratio of the angle of yaw to the angle 
of hank i? always negative. When the airplane is inclined 
to the right (cp positive), it should yaw to the right (j 
negative). Here again, we have the action of the lateral 
component of the weight. 



When is positive a lateral inclination is associ- 

ated with an angular velocity of roll of the same sign tend- 
ing to increase the inclination, the motion then "being un- 
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stable. The ratio %/n 3 is always negative. The ratio 
m 3/^3 reaches rather high absolute values, the angular ve- 
locities of yaw becoming in certain cases considerably 
larger than the angular velocities of roll. 



Bt OSCILLATORY MOTION 
Let p , p , p , p be the maximum amplitudes of the 

1 1 2 ' 1 3 '4 

disturbances of the variables v, p, r , cp associated with 
the lateral oscillation. The ratios of these amplitudes 
are given in the tables below* 



Airplane 


A 


------ — J 


C 


D 


3 


J&L 

vp 4 


1.06 


1 . 345 


1.825 


2.94 


4.28 




1.475 


1 . 45 


1.42 


1.375 


1.347 


% 


1.325 


1.67 


2.545 


3.67 


5.3 


m 


0.9 


1.15 


1.79 


2.67 


3.93 
















— r 



Airplane 



__L. 



i +~ 



2.23 .1.825 j 1.45 
1.715 1.42 ! 1.097 



3.53 I 2.545 I 1..29 



!.075 j 1.79 



1.185 



0.886 



0.622 



0.1785 



0.287 



5 

0.441 
0..329 
0. 20 3 

0.518 



0 . 30 7 



0. 322 



0.1C2 



0.5 65 
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The foregoing tables should "be completed "by the phase 
differences w&i eh determine the sign relations "between the 
disturbances* 



Ai rplane 


A 


B 


C 


D 


E 


»1 ~ ^4 


112° 


109° 


10 3° 


86° 


38° 


^ - ^4 


103° 


104° 


106° 


108° 


109° 


cp — CP 
^3 ^4 


16° 


22° 


15° 


-2.5° 


-50° 


cp - cp 

3 2 


-77° 


« 9° 


-91° 


-110° 

1 


200° 



Airplane 


1 


2 


9 

O 


4 


5 


6 


CP 1 - <P 4 


97° 


10 3° 


108° 


116° 


41° 


56° 


cp a - cp 4 


100° 


106° 


110° 


116° 


82° 


56° 


cp - cp 

3 4 


9° 


15° 


21° 


41° 


240° 


316° 


Cp - Cp 
3 2 


-81° 


-91° 


-79° 


-75° 


+ 158° 


+ 266° 

















The value cp 2 - cp 4 m x indicates that the variable 1 

passes through its maximum of positive amplitudes ?~ T 

seconds "before the oscillation of variable 4. From an ex- 
amination of the tables, the following conclusions nay be 
derived. The disturbance In yaw 6j is, in general , great- 
er than the disturbance in roll 8 cp when the oscillatory 
motion is stable. When it becomes unstable the ratio 
8 J / 8cp decreases. When cp - cp 4 is in the neighborhood of 

90°, the oscillation of the variable 1, that is, of the 
yaw, passes through its maximum positive amplitude about a 
quarter of a period ahead of the variable 4, the lateral 
inclination. The yaw to the left (positive) is followed 
by a bank to the right (positive). 
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The disturbance 8p of the angular velocity of roll 
is not erectly in quadrature with cp. The phase differ- 
ence cp ~ cp 4 varies "between 100 and 110° and is not ex- 
actly equal to 90°. This may he expected since p is not 
exactly the derivative of cp hut 

dep d\[/ 
d = -r-r cos 8 - — cos cp cos G 
dt dt 

The retios P 3 /p 2 result in smaller values than those for 
the spiral mouion. With equal roll the yawing is less in 
the oscillatory motion than it is in the spiral motion. 

Let us now consider the characteristics of the yawing 
motion. When the airplane is stable, c? 3 - cp 4 is included 
between 15° and 40°. The yawing motion passes through its 
positive maximum l/24 to l/9 of a period ahead of the pos- 
itive maximum of cp . The yawing is to the left at the in- 
stant when the left wing is rising, the secondary effect 
of the rotation r thus tending to oppose the roll. The^ 
same conclusion may he drawn from o_ - q> fi « When the oscil- 
lation becomes unstable, for airplanes 5 and 6, the g h '^se 
difference cp - cp- increases; we find 240° and 316 while 
we have cp^ - cp 0 = 153° and 266°. This indicates a revcr- 
sal in the sense of the phenomena since it is easy to see 
that for cp - cp 0 = 180°"" the two rotations r and p aro 

3 ( 4- I V, 

in opposition. At this instant a positive yau r I to tnc 
leftj is associated with a negative roll p (to the right) 
and the secondary effect of the yaw will "be to increase the 
ro 11 • 

The above numerical analysis throws some light on the 
partial motions, of which the resultant motion is composed. 
The characteristics of each of the three motions studied^ 
above depend on the airplane hut are independent of the in- 
i t ial disturbance considered. 

CONCLUSIONS 



1. The analysis of the lateral motion is more compli- 
cated than that of the longitudinal motion. It is thought 
that the results obtained from the complete numerical com- 
putation of examples nay be useful in estimating the effect 
of the different factors. The variable factors here con- 
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sidered were the coefficients of static stability. The mo- 
tions computed were determined for a zero angle of attack. 
It is of course understood that it is not enough that an 
airplane shall "be stable at a low angle of attack. The 
lateral stability should "be investigated, particularly at 
large angles of attack, and this is a problem which at the 
present moment is the most difficult of solution. A log- 
ical continuation of the above study would be the computa- 
tion of the motions for flightt at large angles of attack, 

2 # The theoretical determination of the variables v, 
P» r, qp , and \|/ has already been considered in several 
works. Of these, we may note: 

a) The computations carried out by Halliday on a cer- 

tain airplane (namely, the Bristol "Fighter 1 . 1 } 
for different angles of attack (R. & M. No« 
1306). At low angles of attack, the curves 
found for this airplane have the same general 
appearance as those computed by us for the sta- 
ble airplanes. 

b) An example computed by Mellvill Jones, in volume 7, 

of W. F. Durand's "Aerodynamic Theory." 

REMARKS CONCERNING AUTOMATIC STABILITY 



The present investigation may be considered as prelim- 
inary to the study of automatic stabilizers. We have sought 
to determine first how an airplane of average characteris- 
tics reacts against the principal disturbances it may en- 
counter. It is only after such a study has been undertaken 
that one r.ay inquire as to the effect produced by a stabi- 
lizer on the natural reactions. We have already published 
several ideas on the subject in articles devoted for the 
most part to descriptions of the apparatus proposed. 

Without entering here into the general study of auto- 
matic stabil i zer s , we may point out that the present work 
suggests to us immediately the idea of a stabilizer whose 
sensitive member would be a wind vane or pressure plate. 
The elements here considered as variable were the coeffi- 
cients of static stability - that is, the derivatives of 
the coefficients of the moments with respect to the angles 
of attack and of yaw; these angles may be determined by the 
vane St If the vanes arc utilized so as to operate the cor— 
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trols, an automatic stabilizer is created whose action on 
an airplane will "be analogous to an increase in the static 
stability. 



Let us consider, for 
ity. The moment M is a 
de rivat i vc charac t or i zing 
gravity is: 

dM 
di 



example, the longitudinal stabil- 
function of i and of P. The 
the moment about the center of 



3M 
Si 



aM 

SP 



di 



where dP/di is determined "by the mechanical apparatus con- 
necting the deflection , p of the elevator with position of 
the vane. Our computations are, in fact, referred to the 
total derivative dlvl/di. It scons, from the theoretical 
point of view, all other conditions remaining the same * 
particularly, the total area of the tail surfaces - that 
the production of a moment M "by means of natural static 
stability or by means of an automatic stabiliser would be 
equivalent. The comparison of the curves obtained for dif- 
ferent values of |x shows the effect that may be expected 
from vane stabilisers. The same remark applies to the lat- 
eral stability, whore a vano could be imagined as acting 
on the ailerons or the rtiddcr: 

•dC T SCr dC~ da 
dj 3j 8a dj 

dC^ dCy 3% dY 
dj S Ij + SY dj 

It therefore seems that it should be possible, if an 
airplane possesses certain unfavorable flight characteris- 
tics due to insufficient static stability about one of the 
axes, to .correct this lack of stability by the use of a 
vane stabiliser, thus avoiding a re-design of part of the 
airplane. 



Translation by J. 
National Advi sory 
for Aeronautics. 



Rei s s , 
Commit t oe 
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L01IG-I TUDINAL MOTIONS 
Initial Condition: 6w = -8 and 8© = -0.2 
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LATERAL MOTIONS 
Initial Condition; 6v = +10 m/sec 
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Figure 1.- Scheme of axes. 




Figure 3.- Longitudinal equilibrium. 
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Figs. 5,6 
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Figs, 7,8,9,13,15,16 
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Figure 15.- Components of the disturbance 
in angle of attack. 
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Figs. 10,11 
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Fig. SI 
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Figure 31.- Effeot of an initial disturbance of yaw. 
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Figure 33.- Effect of an initial dieturbance of angular velocity 6. 
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Fig. 26 
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Figure 26.- Distribution of an initial 
disturbance of yaw. 



